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Abstract: 

Working bi-Hamiltonian structure and Jacobi identity in Frenet-Serret frame 
associated to a dynamical system, we proved that all dynamical systems in three 
dimensions possess two compatible Poisson structures. We investigate relations 
between geometry of surfaces defined by potential function of a gradient system 
and its bi-Hamiltonian structure. We show that it is possible to find Hamilto¬ 
nian functions whose gradient flows have geodesic curvature zero on potential 
surfaces. Using this, we conclude that Hamiltonian functions are determined 
by distance functions on potential surfaces. We apply this technique to find 
conserved quantities of three dimensional gradient systems including the Aris¬ 
totelian model of the three-body motion. 

1 Introduction 

A Poisson structure on a manifold is defined by a skew symmetric contravariant 
bilinear form subjected to the Jacobi identity expressed as the vanishing of the 
Scliouten bracket of Poisson tensor with itself IB-0. This structure having no 
non-degeneracy requirement becomes the basic underlying geometry to study 
non-canonical Hamilton’s equations on odd dimensional manifolds as well as the 
Hamiltonian structures of nonlinear evolution equations 0-0- 

The first interesting case of a completely degenerate finite dimensional Hamil¬ 
tonian structure occurs in three dimensions. Many works have been devoted 
to the study of three dimensional dynamical systems with primary concern on 
quantization, construction of conserved quantities, Hamiltonian structures, inte- 
grability problems and their numerical integration using techniques from various 
areas such as Poisson geometry, differential equations, Frobenius integrability 
theorem and theory of foliations 0-|26|. 

In [27] . we reduced the problem of constructing Hamiltonian structures in 
three dimensions to the solutions of a Riccati equation in moving coordinates of 

1 On leave of absence from Department of Mathematics, Yeditepe University, Atasehir. 
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Frenet-Serret frame. All known examples of dynamical systems having two com¬ 
patible and explicit Hamiltonian structures are exhausted by constant solution. 
We concluded that in three dimensions vector fields which are not eigenvectors 
of the curl operator are at least locally bi-Hamiltonian. This structure manifests 
itself in the well-known form of Frenet-Serret triad t = n x b where t is the 
unit tangent vector associated with the given dynamical system and, the normal 
vectors n and b are related to conserved covariants. From the expression of gra¬ 
dient operator in Frenet-Serret frame (c.f. Eq.([3j)) we observe that unit tangent 
vector arises from the Cartesian gradient of a function of arclenght variable. 

In this work, we shall continue to investigate the local structure of bi- 
Hamiltonian systems in three dimensions assuming that they are described by 
gradient vector fields. We shall relate ingredients of bi-Hamiltonian structure 
to geometry of surfaces described by the potential function F. In particular, 
we shall prove that there are two Hamiltonian functions related with geodesic 
coordinates on the potential surfaces. 

1.1 Content of the work 

In the next section, we shall first develop differential calculus in Frenet-Serret 
frame that will be used throughout the paper. Then, we shall survey on con¬ 
ditions for construction of Frenet-Serret frame for a given vector held in R 3 to 
ensure its existence. This will extend our previous result in m that excludes 
eigenvectors of curl operator. 

In Section 3, we shall summarize ingredients of bi-Hamiltonian systems in 
three dimensions. In particular, we shall identify a Poisson bi-vector with a 
locally integrable (in the sense of Frobenius) vector held in three dimensions 
and, will work with the latter. We shall first reduce the Jacobi identity into 
Riccati equation and then, assuming we have independent solutions, we shall 
exhibit relations between conserved Hamiltonians and Poisson vectors. 

In Section 4, we shall start with geometric characterization of potential sur¬ 
faces in the normal coordinates of the Frenet-Serret frame. By considering gra¬ 
dient hows of Hamiltonian functions on potential surfaces, we shall prove that it 
is possible to find Hamiltonian functions whose gradient hows on the potential 
surface have geodesic curvature zero. Using this result, we shall conclude that 
Hamiltonian functions are determined by distance functions on potential sur¬ 
faces. Finally, we shall prove that Hamiltonian functions of the gradient system 
are related with geodesic coordinates of the potential surfaces. 

In Sections 5, we shall present examples of gradient dynamical systems which 
are bi-Hamiltonian and, work out in details the geometry of potential surface. 


2 Frenet-Serret Frame 

Let t (x,y,z) be a given unit vector held in R 3 endowed with Cartesian coor¬ 
dinates x = ( x , y 1 z). We may assume t to be a unit tangent vector to a curve 
t —> x(f) in R 3 . This may be the solution of dynamical system in Ea. (fT51) . In 
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this case, we can choose t = v/||v||. Locally, one can always lift t to an or¬ 
thonormal frame in R 3 in infinitely many ways. Let (t, n, b) be such an arbitrary 
orthonormal frame satisfying 

t = n x b, n = bxt, b = t x n. 


We introduce the directional derivatives along the triad (t,n, b) as 


d s = t ■ V d n = n ■ V d b = b ■ V (1) 


so that the variables ( s,n,b ) are the coordinates associated with the Frenet- 
Serret frame. Assuming the Cartesian coordinates are functions x = x(s,n.b) 
of Frenet-Serret coordinates we find, using Eq. © the Jacobian matrix 


d(x,y,z) 
d(s, n, b) 


d s x 

d n x 

d b x 

d s y 

9nV 

d b y 

d a z 

d n z 

d b z 


t n b 


whose determinant 


det 


( d(x,y,z) \ 
\ d(s, n, b) ) 


t(n x b) = 1 


is non-zero and hence the inverse transformation 


s = s(x), n = n(x), b = 6(x) (2) 

exists locally, that is, in a sufficiently small neighborhood of a given point xq £ 
R 3 . These functions may be obtained by integrating the quantities 


ds = t • dx, dn = n • dx, db = b ■ dx 


the last two of which implies n =constant and b =constant when restricted to 
the curve x(f). 

2.1 Differential calculus 

By inverting equations © we get the expression 

V = td s + n d n + b d b (3) 

for the Cartesian gradient in Frenet-Serret frame. For future reference, we define 
the helicities m 


Ht = t • V x t, Hn = n • V x n, Hb = b • V x b (4) 

and the cross-helicities 

Utn = t ■ V x n, H nt = n ■ V x t, U nb = n ■ V x b (5) 

Htb = t • V x b, Hbt = b • V x t, Hbn = b • V x n. 
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which measure holonomicity of Frenet-Serret triad. From the coefficients of the 
basis vectors in expansions of curls into (t,n, b) we obtain 

V x t = H t t + Hutu + H bt b 

V x n = HtnX + H n n + H bn b (6) 

V x b = Htb t + 'Hnbn + Hbb 

and from the orthonormality of basis vectors we have the divergences 

V • t = Hbn — H n b 

V-n = Htb-Hbt (7) 

V • b = Hnt. ~ Ht.n- 

Using the vector identity 

V(u ■ v) = (u ■ V)v + (v ■ V)u + u x (V x v) + v x (V x u) (8) 
we compute the derivative of t in direction of t to be 

d s t= (t • V)t = —t x (V x t) = — (n x b) x (V x t) = n Hu — b H n t- 
Repeating for the other directions we have the derivatives 

d s t = nHbt ~ b Hnt 

d n n = b Htn - W-bn (9) 

d b b = t 'Hub ~ nHtb 

of basis vectors along their respective directions. To find other derivatives, add 
the identity 

V x (u x v) = (v • V)u — (u • V)v + (V • v)u — (V ■ u)v 
to Eq.© to obtain 

2(v • V)u = V(u • v) + V x (u x v) — (V • v)u + (V • u)v 

— u x (V x v) — v x (V x u) 

and let u = t, v = n to have 
2(n ■ V)t = 2d n t 

= V x b — (V ■ n)t + (V ■ t)n — t x (V x n) — n x (V x t). 

Last two terms may be expressed as 

t x (V x n) = (n x b) x (V x n) 

= H n b - H bn n 

n x (V x t) = (b x t) x (V x t) 

= Hbtt — Htb. 
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Collecting these for d n t and repeating similar computations for other derivatives 
we obtain 


d n t = 

Tibn™ + -('Ht — Tin + Tib) b 

d h t = 

--(Tit + Ti n — Ti b )n — Tinbb 

d t n = 

-Tibtt + \{Tit-Ti n -Ti b )h 

d b n = 

-("Hj + Tin — Tib) t + Titbb 

d t b = 

Tintt — -( Tit — Tin — TTb)m 

On b = 

——{Tit — Tin + Ti b ) t + Titn n. 


2.2 Constructing Frenet-Serret frame 

The choice of orthonormal frame is determined by the choice of n. In [57], we 
introduced such a frame assuming that the unit tangent is not an eigenvector 
of the curl operator. Here, we want to release this assumption and prove the 
existence of an orthonormal frame for all smooth dynamical systems in three 
dimensions. Our result will rely on the eigenvalue problem 

V x t = A(x)t (11) 

for the curl operator. If t is not an eigenvector of the curl operator we have 
(V x t) x t / 0 and we recover the result of [[27]. If however, the eigenvalue 
equation m holds, then Tit = t • V x t =A(x). At each point x, the eigenvalue 
A(x) will define a surface with normal V"Ht(x) if ’Ht(x) is not a constant function. 
We distinguish two cases depending on whether the unit tangent has components 
lying on this eigensurface or not. If it has, then we choose the normal on the 
eigensurface. If t is completely aligned with the surface normal, then we recall 
a result of Chandrasekhar and Kendall in [28] that there exist a constant unit 
vector defining an eigenvector of the curl operator and construct a frame with 
this constant unit vector. In the remaining case with A(x) = 0, we have a 
surface whose gradient is the unit tangent and we choose the frame using lines 
of curvature of this surface. More precisely, following result proves that there 
are canonical liftings to Frenet-Serret frames. 

Proposition 1 Given a nonzero vector field v £ R 3 , let t = v/||v||. Then, the 
vector field n can be chosen as follows 
1. // (V x t) x t ^ 0 then let 


n = 


(V x t) x t 
(V x t) x t 


and we have necessarily Tint = n ■ V x t = 0. 


( 12 ) 
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2. If (V x t)xt = 0, then V x t=77*t and we have necessarily H n t = Hbt. = 
0, or equivalently d s t = 0. We distinguish two cases: 

2a. if VW( xt^O then choose 


VH t x t 

WHt xt 


(13) 


Ufri. */ V"H f x t = 0 and Ht 
unit vector a swc/i that 


= constant ^ 0 , t/ien f/iere exists a constant 


a x t 

n = -- 

a x t 


(14) 


and hence H n t = Hbt = 0. 

2bii. if S7Ht x t = 0 and Ht = 0, then there exists a surface with normal 
t and the Frenet-Serret frame is the Darboux frame of the lines of curvature on 
this surface. We have Hnt = Hbt — 0. 


Proof. Case 1. follows easily. For case 2., assume (V x t) x t = 0. Then, V x t 
is proportional to t, V x t=Htt and we have H n t = Hbt = 0 which are the 
coefficients of unit vectors in the expression for d s t. In order to construct the 
normal vector n we have two subcases depending on whether S7Ht is non-zero 
and parallel to t. If X7Ht xt/0 then define the unit normal as in Ea. dTdl) . If 
V"Ht x t = 0 we have VHt proportional to t. If the proportionality function is 
zero, we have Ht =constant. In this case, t is an eigenvector of the curl operator 
with constant eigenvalue Ht- By a result obtained in [28] there exists a scalar 
function %j) satisfying 

Atp + Hfip = 0 ( 15 ) 

and a constant vector a such that the eigenvector t can be expressed as 

t = — V x fipSL + V x (V’a)) • (16) 

Ht 

Then, we define the normal vector by Eq. (flTl) . It follows that 

H nt = n • V x t = 77———777 ‘ Htt = 0 
11 a x 111 

Hbt = b ■ V x t = (t x n) ■ Htt = 0. 


If the proportionality function is non-zero, taking curl and then dot product 
with t we get Ht = 0 which is the integrability condition for the unit tangent 
vector. Since, this is a subcase of 2 with V x t = HtX 1 we have V x t = 0. 
Locally, there exists a function F (x) such that t = VF (x). Choose n to be the 
unit tangent vector of the line of curvature of this surface, namely, a vector n 
satisfying 

n-(t x (n-V)t) = 0 (17) 

which implies, from Eq. lfTUl) Ht = 0, H n = Hb■ Moreover, we have H n t = Hbt = 
0 as in the previous case. ■ 
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3 Bi-Hamiltonian Structure 


We shall summarize the necessary ingredients of the bi-Hamiltonian formalism 
in three dimensions. See 0-126] for details and examples. For x = {a; 1 } = 
(. x,y,z ) £ I 3 , i £ I and overdot denoting the derivative with respect to t, we 
consider the system of autonomous differential equations 


dx 

dt 


= v(x) 


(18) 


associated with a three-dimensional smooth vector field v. Eci. (fl8l) is said to 
be Hamiltonian if the right hand side can be written as v (x) = H (x) (dH (x)) 
where H (x) is the Hamiltonian function and f l (x) is the Poisson bi-vector 
(i.e. a skew-symmetric, contravariant two-tensor) subjected to the Jacobi iden¬ 
tity [fl(x),fl(x)] = 0 defined by the Schouten bracket pQ. In coordinates, if 
di = d/dx l , the Poisson bi-vector is f l (x) = ( x) dj A dk , with summation 

over repeated indices, and the Jacobi identity reads = 0 where [jkl] 

denotes the antisymmetrization over three indices. It follows that in three di¬ 
mensions the Jacobi identity is a single scalar equation. One can exploit vector 
calculus and differential forms in three dimensions to have a more transparent 
understanding of Hamilton's equations as well as the Jacobi identity. Using the 
isomorphism 

Ji = e ijk W k i,j,k = 1,2,3 (19) 

between skew-symmetric matrices and (pseudo)-vectors defined by the com¬ 
pletely antisymmetric Levi-Civita tensor £ijk, we can write the Hamilton’s equa¬ 
tions and the Jacobi identity as 


v = J x Vff J • (V x J) = 0, 


( 20 ) 


respectively. In this form the Jacobi identity is equivalent to the Frobenius 
integrability condition J Ad J = 0 for the one form J = Jidx 1 . It is the condition 
for J to define a foliation of codimension one in three dimensional space inra- 
m- a distinguished property of Poisson structures in three dimensions is the 
invariance of the Jacobi identity under the multiplication of Poisson vector J (x) 
by an arbitrary but non-zero factor. The identities 


J • v = 0, Vif • v = 0 (21) 

follows directly from the Hamilton’s equations in (1201) . The second equation 
in m is the expression for the conservation of Hamiltonian function. A three 
dimensional vector field v (x) is said to be bi-Hamiltonian if there exist two 
different compatible Hamiltonian structures mm- In the notation of equation 
(Ullh this implies 

v = Ji x VH 2 = J 2 x VHi (22) 

for the dynamical equations. The compatibility condition for Ji and J 2 is 
defined by the Jacobi identity for the Poisson pencil Ji + CJ 2 for arbitrary 
constant c. 
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3.1 Jacobi identity in Frenet-Serret frame 

Given a vector field v, v/ ||v|| is the unit tangent vector t to the flow of v. It 
follows from the identity J • v = 0 that the Poisson vector J has no component 
along the unit tangent vector t. Hence, we set 

J = An + Bb (23) 

for unknown functions A (x) and B (x) satisfying A 2 +B 2 ^ 0. Assuming A ^ 0 
and defining the function p = B/A the Jacobi identity for J =A(n + //b) reduces 
to the Riccati equation 


dsfX = 'Hn + n{T~Lnb + Bbn) + p 2 TLb (24) 

and 

d a ln A = d a \n.\\v\\- (Mb-Hub (25) 

in the arclenght variable s. The Riccati equation (1241) is equivalent to a linear 
second order equation and hence possesses two linearly independent solutions 
leading to two Poisson vectors for dynamical system under consideration. The 
Hamiltonian form of dynamical equations implies that the Poisson vectors ob¬ 
tained from solutions of Riccati equation are always compatible. Thus, we 
conclude that 

Proposition 2 All dynamical systems in three dimensions possess two compat¬ 
ible Poisson vectors. 


3.2 Poisson vectors and conserved quantities 

Once we have the independent solutions p 1 and p 2 °f the Riccati equation, we 
can form the compatible Poisson vectors 

Ji=Ai(n + /qb), J 2 =A 2 (n + p 2 b) (26) 

with conformal factors A\ and A 2 . The construction of corresponding Hamil¬ 
tonian functions to form a bi-Hamiltonian pair requires integration of these 
Poisson vectors. 


Proposition 3 The conserved covariants for Ji and J 2 are 


VHi 


A 1 A 2 {p 2 - /q) 


Ji, ViJ 2 = 


A 1 A 2 {p, 2 - Mi) 


respectively. 


(27) 


Proof. By definition, the conserved Hamiltonians satisfy v-Vifi = v-Vi/i = 0. 
That means, the gradients lie on the space spanned by {ni, n 2 } or, equivalently, 
by { Ji, J 2 }. We, thus, form the linear combinations 


W — aJi~\-bJ 2l V H 2 — cJ\+dJ 2 
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and determine the coefficients. Bi-Hamiltonian form in Ea. tB^l) together with 
the orthonormality of the basis {t, n, b} imply 


A 1 A 2 {ii 2 - Mi)’ 

The identities VxViJi = Vx VH 2 = 0 dotted with Ji and J 2 result, recalling 
also the compatibility condition, in the fact that a and d are also multiples of 
this function. It turns out that, for each Hamiltonian, there is an equivalence 
class of Hamiltonian functions whose gradients differ by the functions a and d. 
Hence, without restriction to generality, the conserved covariants are given by 
Eq.©. ■ 


4 Gradient Systems 


We assume that there exists a potential function F for the velocity field 

dx , _„ , . , 

— = v(x)=VF(x) (28) 

of a given dynamical system in three dimensions. In this case, we have V x 
v (x) = 0, so that 

V x t = t x V In ||v|| (29) 

and hence, t • V x t = 0. Then, either (V x t) x t / 0 or V x t = 0. In the 
first case we take n = ((V x t) x t)/||(V x t) x t|| with n • V x t = 0 as well. 
Thus, the Frenet-Serret frame for a gradient system with potential function F 
may consists of unit vector fields 


v VF (Vxt)xt u Vxt 

vf[ ~ 11V W| | ’ n ~ || V x 111 ’ ~ || V x t 


(30) 


In the second case, there exists a surface with normal t, which is nothing but any 
potential surface of F (x) and the Frenet-Serret frame is the Darboux frame of 
lines of curvature on the potential surface. In other words, t will be unit normal 
of the potential surfaces and choosing the Frenet-Serret frame as the Darboux 
frame of the lines of curvature on the potential surface will work in both cases. 

Now, since {t, n, b} forms an orthonormal frame, the arclengths of their 
integral curves (s,n, b) form a local coordinate system around any point p € 
K 3 provided that v (p) 7 ^ 0. On the other hand, 


implies that 



v t = VF (x) 


(31) 

t-VF (x) 

= d s F (s, n, b) = 

IMI 


n-VF(x) 

= d n F (s, n, b ) = 

0 

(32) 

b-VF(x) 

= d b F(s,n,b ) = 

0 
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Hence F (x) = F ( s ) in the coordinate system defined above. Therefore, with 
the assumption that v (p) = X7F (p) ^ 0 the potential surfaces are 

F (x) = F (s) = d (33) 

which implies s = c. As s is the arclength of integral curve of E51) . its flow iv 
acts on the s coordinate and therefore on the potential surfaces by translation, 
i.e. $ 0 - (s, n, b) = (s + a, n, b). If H (x) = H (s, n, b) is a Hamiltonian function 
for (l28l) . then it will be invariant by this flow H (s + a, n,b) = H (s, n, b) , which 
implies that H is independent of s, in other words fixing s = c and letting a = s 
lead to 

H (s,n, b) = H(s,n, b) | s=c . (34) 

The invariance of Hamiltonians under the flow implies that it is sufficient to 
determine the Hamiltonian function on any potential surface. Namely, a Hamil¬ 
tonian function can be reduced to a Hamiltonian function on potential surfaces. 
However, the converse may not be so straightforward. For the construction of a 
Hamiltonian function out of a function defined on potential surfaces, or equiv¬ 
alently, to extend a gradient vector X/H on the potential surface to a gradient 
vector on R 3 , the vectors t and t x Vf?, which are perpendicular to Vi?, must 
be tangent to a surface in R 3 . This condition can be written as 

[crt-V, (/?t x Vi?) • V] = 0 (35) 

for some functions a and /3 on R 3 , where [•,•] denotes the bracket of vector 
fields. Our purpose is to prove that 

Theorem 4 The Hamiltonian functions of a gradient system defined by a po¬ 
tential function F are determined by geodesic distance functions defined by non¬ 
conjugate points on the potential surfaces. 

To obtain this result, first we are going to show that, geodesic distances on a 
potential surface are defined by a gradient system. Then, we will prove that all 
Hamiltonian functions on potential surface are generated by geodesic distances 
to non-conjugate points. Finally, we will show that these distance functions can 
be extended to Hamiltonian functions on R 3 . 

4.1 Differential geometry of potential surfaces 

We give geometric parameters of potential surface and of an arbitrary curve on 
it. 

Proposition 5 The fundamental forms of surfaces F (s) = c are 
ds\ = dn 2 + db 2 

ds\ = Hbndn 2 + ('Hb—'H n )dbdn —'Hnbdb 2 . 
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IfX (n (a ), b ( a )) is a curve, parametrized with the arclenght a, on the potential 
surface, then the normal and the geodesic curvatures are 

n n = t; 2 'Hbn+t;y(7-Ln—'Hb) + V~'Hnb (36) 

Kg = db€ - d n y - i'Htn-y'Htb, (37) 

where we denote 

dn/da = idb/da = y. 

Proof. Let X = X(c, n, b) be a potential surface of the function F determined 
by a constant c. The tangent vectors and the unit normal vector of this potential 
surface are 

<9„X = n, d b X = b, a„Xx d b X = t. 

The first fundamental form is obviously ds\ = dn 2 + db 2 . To find the second 
fundamental form, first note that the normal vector t(c, n, b) of X(c, n, b) is 
independent of s. We compute 

dt-dX = (d n tdn + d b tdb) ■ (n dn + bri&) 

= n • {d n tdn 2 + d b tdbdn) + b • (d n tdndb + d b tdb 2 ) 

= (n • d b t + b • d n t)dbdn + n • d n tdn 2 + b • d b tdb 2 
= 'H bn dn 2 + (' H b —'H n )dbdn — 'H nb db 2 


and thereby obtaining the second fundamental form. For curvatures, the unit 
tangent vector of the curve X (n (a ), b (a)) is 


_ dn db, 

T = ^ = ^ n+ ^ b 


dX 

da 


with ( dn/da ) 2 + ( db/da ) 2 = 1. We compute the curvature vector 

clT df .dn dn db 

— = -f-n+£— + -r^b+r?— = kN 
da da da da da 

= (£d n £ + yd b £)n + £(£d„n + yd b n) 

+ {£d n y + yd b y)b + y(£d n b + yd b b). 


(38) 


The second and the third parenthesis can be expressed, using formulas from 
calculus in Frenet-Serret frame, as 

fd n n + pd b n = £(H tn b-'Hbnt) + y{-U t n^+-{'H n -'Hb)t), 

fd n b + rid b b = y(-'H tb n+'H n bt) + ^('H tb i>+-{'H n -'H b )t), 

with which the curvature vector becomes 

dT 

—r~ — (—£ 2< H6n+£T7CH n — Mb) + V 2 ^nb) t 

da 

+(-yd n y + yd b £ - fy'Htn-tf'H.tb )n 
+(£,d n y - fd b £, + f 2 'Htn+fy'H tb )b 

= {-fn bn +^{n n -n b ) + y 2 u nb ) t 

+ (<9&£ - d n y - Z'Htn-y'Htb)(yn-fb) 
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where we used = — qd n p and qdbP = —fdbf. This gives the decomposition 

dT , 

— = K n t + n g {i] n-£b) 

of the curvature into tangential and geodesic components. ■ 

A curve parametrized with arclength is a geodesic on a surface if n g = 0 
for all points on the curve. We will see that this condition is necessary for the 
integrability of forms forming transformations to geodesic coordinates on the 
potential surface of a dynamical system. There will be two such coordinate 
systems, each corresponds to construction of a Hamiltonian function on the 
surface. 


4.2 Flows of Hamiltonian functions on potential surfaces 
and geodesic distances 

In this section, we shall prove two propositions which relate the flows of Hamil¬ 
tonian functions on potential surfaces with geodesic distances. First proposition 
tells that geodesics are integral curves of gradient flows of geodesic distances on 
potential surfaces. Second proposition proves that geodesic distances defined 
by non-conjugate points, i.e. points which are joined by a unique geodesic, 
are functionally independent and therefore defines two Hamiltonian functions. 
Hence these functions generates all Hamiltonian functions that can be defined 
on the potential surface. 


Proposition 6 Finding two Hamiltonian functions for the gradient system VP 
on the potential surface F = c amounts to finding geodesic distances on the 
potential surface. 


Proof. Let (it, v) be orthogonal coordinates on the potential surface F (x) = c 
for a given Riemannian metric 


(9ij (u,v)) = 


9uu (u, v) 0 

0 g vv (it, v) 


and let (q,p) be another orthogonal coordinate system such that 


(Gij (<1,P)) = 


Gqq (pp) 

0 


0 

G pp ( q,p ) 


Above two metrics are related by 


9uu v) — Gqq (<?, p) Q u + Gpp (q, p) p u 

g vv (it, u) = Gqq ( q , p) q 2 v + G pp (q, p) p 2 v (39) 

0 = G qq (q, p) q u q v + G pp ( q , p) p u p v 

Note that, although ( u,v) and (q,p) are orthogonal coordinates on the surface, 

the coordinate transformation between them need not be orthogonal, that is, 
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the Jacobian determinant 


J = QuPv - qvPu = 

is not necessarily equal to unity. Defining 


QuuQvv 

GqqGpp 


= 
/ u 


ll = 


99 


Qvv 

G~ 


Tu = 


99 


G 


pp 


Iv = 


Qvv 

G~ 


(40) 


(41) 


pp 


and introducing a parameter 9, we can write equations dMD as a system of 
differential equations 


dq 

du 

dq 

dv 


ll cos 9 , 

-7®sin0, 


dp 

du 

dp 

dv 


= lu sin 9 


= 7 ^ cos 9. 


(42) 

(43) 


for transformations between two orthogonal coordinates. Equations (1421) and 
(l43l) are subjected to two integrability conditions 


q . „dd 9 a 99 

7^0^- 7 ; cos 0 ^ 

„ „ae „ . „ae 

TScose^+TSsmS^ 


d^ q u d^ q . 

1 - COS 9 + —— sm 9 

dv du 

&S cosS _&i sin «, 

ou ov 


(44) 

(45) 


The characteristic curve of the first integrability condition in Eg. (PHI) is the 
integral curve of the dynamical system defined by 


du 

dt 

dv 

dt 

d9 

dt 


—"/l cos 9 


ll sin# 
dll 


dv 


cos 9 


dll 

du 


sin 9 


(46) 


and the arclength a\ of this integral curve satisfies 


da i 
dt 

In the transformed coordinates ( q,p ), the first two equations for the character¬ 
istic curve become 



dq 1 da i dp 

dt \jG qq dt 1 dt 


( 47 ) 


hence, the arclength p ± in the transformed coordinates satisfies 


dp-y da i 
dt dt 
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and Eq. (1471) in arclength parametrization becomes 


-T(! ) = -v®7 


dp 1 \ P 


The geodesic curvature of the characteristic curve q = c is given by 

1 9 In G n 


*qq 


2^/G~ dq 


( 48 ) 


(49) 


Similarly, for the second integrability condition (1451) , the characteristic equa¬ 
tion is the integral curve of the dynamical system 


du 

dt 

dv 

dt 

d6 

dt 


— = lv sin 9 

lu cos 6 
djP 


(50) 


. cos 9 — sin 9 
ou ov 


and the arclength ct 2 of this curve is determined by 

da 2 


dt 


— \J flW'Yu- 


In the transformed coordinates ( q,p ), the first two equations for characteristic 
curve are 

dq n dp 1 da 2 
dt 


= 0, 


dt dt y/G P p dt 

and the arclength p 2 in the transformed coordinates satisfies 


(51) 


dp 2 

dt 


da 2 

dt 


Equation (1511) in this parametrization becomes 

4r( ! ) = -Vg^ 


dp 


P 


The geodesic curvature of the characteristic curve p = c is 

1 d In G vv 


^\fGqq d P 


(52) 


(53) 


(54) 


To arrive at the conclusion of the proposition we will restrict the character¬ 
istic curves to be geodesics. To this end, we note that we can choose two types 
of geodesic coordinates on a surface. One transforming the metric as 


3uu (u, v) 0 

0 g vv (u, v ) 


1 0 
0 Gpp (q,p) 


(55) 
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whereas, the other transforming the metric as 

( 9uu{u,v) 0 \ _ / G qq (q,p) 0 \ 

V 0 9w (u, v) ) \ 0 1 ) ' 

The first transformation reduces the equations fl51) and (HU to 



for the characteristic curve of the first integrability condition in Eq. (l44l) and its 
geodesic curvature. Similarly, the second transformation reduces the equations 
(|5dl) and (|54|) to 



which are associated to characteristic curve of the second integrability condition 
in Ea. (l45l) . 

From Ea. (l47l) and Eq. m we have 

dq da i dp da -2 

dt dt ’ dt dt 

and without restriction of generality, we may assume that 

Hi = —q = ai H 2 = p = a 2 (59) 

which are arclengths along geodesic curves, namely geodesic distances. ■ 

Next, we are going to prove that geodesic distance functions defined by non¬ 
conjugate points are two Hamiltonian functions on the potential surfaces. 

Proposition 7 The Hamiltonian functions of a gradient system defined by a 
potential function F are determined by geodesic distance functions on potential 
surfaces. 

Proof. Since level surfaces, F (x) = c, are closed subspaces of R 3 , which is a 
complete metric space, level surfaces with the induced metric are also complete. 
Then, by Hopf-Rinow theorem m they are also geodesically complete, which 
implies the existence of a geodesic between any two points on a level surface. 
Choose and fix a level surface F (x) = c, which we denote by S c , and a point 
Pi £ S c - Then the first Hamiltonian function on the potential surface, Hi (x), 
can be constructed as 

Hi (x) = d (x, pi) = The length of the geodesic joining x GS C to pi. 

Since the gradients of distance functions have unit norm, 

||Vi7i (x)|| = ||Vd(x, pi)|| = 1 for x GS C 

and belong to tangent plane of S c at x, it is a Hamiltonian function of (l28l) . 
whose gradient has unit norm on the potential surface. 
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For the second Hamiltonian, let T (pi) be the set of intersection of all 
geodesics through pi. Choosing another point p 2 ^ T on the level surface S c , 
and repeating the same construction, the second Hamiltonian function is ob¬ 
tained. This second function will be independent from the first one. In the 
rest of the proof we will show this. Assume, on the contrary, that the second 
distance function is not independent. Then, 

Vd (x, pi) x Vd (x, p 2 ) = 0 


which implies that 

Vd (x, p 2 ) = AVd (x, pi) 
and since they have unit norm 

II Vd (x, pi) || = || Vd (x, p 2 ) || = 1 
we obtain A = ±1. Now, if A = 1 then, 

Vd (x, pi) = Vd (x, P 2 ) 


and therefore 

d(x, p 2 ) = d(x, pi) +k. 

Choosing x as the midpoint of the geodesic combining pi and p 2 implies that 
k = 0 on the surface, and therefore p 2 = Pi 6 T (pi) which contradicts our 
assumption. On the other hand if A = — 1, then 

Vd (x, pi) = — Vd (x, p 2 ) 


and therefore 


d(x, pi) + d(x, p 2 ) = k 

Setting x = pi implies that k = d (pi, p 2 ), and the condition becomes 
d(x, pi)+ d(x,p 2 ) =d(pi,p 2 ). 


By triangle inequality, this means that the points x, pi,p 2 lie on the same 
geodesic for all x on the level surface. Hence, Vd(x, pi) = —Vd(x, p 2 ) is 
possible only if every geodesic through pi contains p 2 , and therefore p 2 G T (pi) 
which again contradicts to our assumption. Namely, distance functions defined 
by these two points specified above are functionally independent. Since the 
gradient of these two functions span the tangent plane at any point, gradient 
of any third function obtained in this way will be linearly dependent on the 
previous two, therefore it will be a function of H\ and id 2 . ■ 

To sum up, if we take X(n(er,:) ,&(oy)) to be the gradient flow for i — th 
Hamiltonian function on the potential surface, we get 


£ 


da 


db 0 u 
r t = ~r = VbHi, 
da 


i = 1,2. 
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together with 


\\S7Hif = (d n Hi) 2 + (dbHi) 2 = 1 

Furthermore, from Eas. (l26l) and (l27l) we can identify the partial derivatives of 
Hamiltonian function and obtain 

d b Hi = n,d n H t (60) 

or rj = and rj 2 + £ 2 = 1. Then, the gradient of Hamiltonian functions on 
potential surface become two different representations of the unit tangent vector 
of X (n (a) , b ( a )) 

T = X7H l = ^=L=(n + f , l h) = ^, < = 1,2. (61) 

y/l + tf ll J dl 

The integrability condition 

dn {dbHi) - db {d n H t ) + UtndnH + HtbdbH = 0 (62) 

of Eq.® 

Kg = --- 3 - (dnH + ,n + (l + /r 2 ) C Htn + iMtb)) = 0 (63) 

(1 + M 2 ) 5 

is the vanishing geodesic curvature described in Eo. (l37l) . 

Note that the Riccati equation (l24ll describing /i (s, n, b) determines the par¬ 
tial derivative with respect to s variable and allows an arbitrary dependence on 
n and b variables, while the vanishing geodesic curvature condition depends only 
on the derivatives with respect to n and b. Therefore, the choice of geodesic dis¬ 
tances as Hamiltonian functions puts a restriction on the arbitrariness admitted 
by Eg.®. 

4.3 Local extension of Hamiltonian functions on potential 
surfaces 

On any potential surface, F (x) = c, we had two gradient vectors Vffi, VJJ 2 and 
the unit normal of the surface t. One can extend these Hamiltonian functions 
on potential surfaces to functions on M 3 in infinitely many ways. Geometrically 
this amounts to embedding geodesic curves to surfaces in K 3 , which are the 
potential surfaces of Hamiltonian functions of dynamical system (12811 under 
consideration. These surfaces must contain integral curves of Ea. ([28]) . In other 
words, tangent plane of surface obtained by extension of Vi?; must contain all 
vectors perpendicular to each Vifj in M 3 , in particular, t and Vi4j x t. In fact, 
since t and S/Hi x t are linearly independent by definition, they span a two 
dimensional plane. The following proposition proves that they always integrate 
to potential surface of the corresponding Hamiltonian function. For convenience 
we will assume that Ea. (l28l) is of the form 

v (x) = VF (x) = Ji (x) x S7H-2 (x) = J 2 (x) x VTJi (x) (64) 
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of a bi-Hamiltonian system, where 


3 1 =4>VH 1 , 3 2 = -<jNH 2 


(65) 


and the conformal factor is given by 

0 _ A 1 A 2 (^ 2 — /r 1 ) 

IMI 


Proposition 8 


( 66 ) 


t-V (Jj x t) • V 


= 0, 


* = 1,2 


(67) 


Proof. Let w = £n + 77 b be a vector field in the tangent plane of a potential 
surface of the function F. Then 


[t-V, w-V] = -?4t£t-V+ (d s £ - Hbrit, + 'Hn'n) n 'V 
+ ( d s r] - H b £, + HubV) bV. 


( 68 ) 


If w = Vi? for some function H on the potential surface, we get 

[t-V.VLT-V] = -'H b td n m-\7-(2'HbndnH + (H b -'H n )d b H)n-\7 

+ (2H nb d b H+{H n -H b )d n H)b-V. 1 j 

Using (l69l) . it is easy to compute 

[t-V, (Vif x t) -V] = -HbtdbHt-V- (V • t) (VET x t) -V (70) 


which is sufficient to show that the space spanned by t and VIE x t is tangent 
to a surface in R 3 . Using the equations 

dgiii = Wn + n i {‘H nb + V-bn) + i — 1,2 (71) 

d s \aAi = <9 s ln||v|| - ^{H b - U nb * = 1,2 (72) 

defining the Poisson structures we obtain 

ds 111 (/Z 2 - Ml) = (Hnb + ?tfm) + (a* 2 + Ml) % (73) 

In py = - 'Hnb- (74) 

Eas. (173l) and (17-El) when used in Eg. (IfETH) result in 

8 S In = Ubn - U nb = V • t (75) 

IMI 

for the divergence of t. This leads, with Eq. dTUl) . to 


t-V 


(Vi/i x t) -V 


-^--UbtdbH, t-V 

V 


(76) 
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(77) 


By Eg. (1711) we have 
0 n ifi = - —, 


d b H 1 = -^, d n H 2 = % d h H t = ^ 


and we can rewrite (17611 in the following form: 


A, 


t ' V ^Ti — iT (Mi n - b ) -V 




"Hbtt-V 


Now, for any vector field u satisfying 

[t-V, u-V] = At-V 

we can find an integrating factor a such that 

[at-V, u-V] = (a\ — u-Vq) t-V= 0 

and therefore 


Taking 

Eq.© amounts to 


aA — u-Va = 0 


u = 77—77 (/^n — b ) and A = 

v v 


d b a - Hid n a = / 1 / Hbt * = 1 , 2 . 
Since /i 2 — Pi yf 0 we simply get 

dp.OZ — 'Hbti d b CX — 0 . 


(78) 

(79) 

(80) 
(81) 
(82) 

(83) 

(84) 


To find a solution to (1511) . note that, for any function ip (s ), the derivative 
d s ip(s) satisfies 


d n d s ip(s) = HbtdsPis) 
d b d s p (s) = -Untdsp (s). 


(85) 

( 86 ) 


According to Proposition 1 above, % nt = 0 for the cases 1 and 2bii of Frenet- 
Serret frames that can be applied to gradient systems. Therefore, we have 


and choosing 


d n d s ip (s) = HbtdsP (s) 

d b ds<p (s) = 0 


1 


(87) 


( 88 ) 


d s p (s) 

solves Eg. (1511) . It is possible to get a more specific solution using Eos. (l32l) and 
m which state that 

||v||=0 s F(s). (89) 
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Therefore one may choose 


a = 


1 

^1 


(90) 


which proves the proposition ■ 

Above proposition proves that the arclengths of integral curves of vector 
fields 



provide parametrizations and hence a local coordinate system for the potential 
surface Hi = c. 


5 Examples 


5.1 Sphere 

Let r = ( x,y,z ) and r = ||r|| = x 1 + y 2 + z 2 . Consider the gradient system 


dr (t) 
dt 


r (*) 


(91) 


with the potential function 


F( r) = 


(92) 


where level surfaces are spheres. To construct the Frenet-Serret frame, we begin 
with the unit tangent vector 



(93) 


and take normal and binormal vectors as the lines of curvature of level surfaces 
of the potential function. Since level surfaces are spheres, which are surfaces 
of revolution, their lines of curvatures are latitudes and longitudes. Adapting 
spherical coordinates 

x = r sin 0 cos 9 

y = rsin0sin0 (94) 

z = r cos 0 


the Frenet-Serret frame can be written as 


t = sin 0 cos 9\ + sin 0 sin 6 j + cos 0k 
n = — sindi + cosdj (95) 

b = cos0cos0i + cos0sin6*j — sin 0k. 

For the construction of Hamiltonian functions, we will fix the unit sphere. The 
distance function on the unit sphere is 


d (. Pi,P ) = arccos (P,-P) 


(96) 
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where P and P; are position vectors of the points P and Pi. Written explicitly, 
if 

P* = (sin^cos^sun^sinf^cos^) (97) 

Q = (sin 0 cos 0, sin</>sin0,cos<?!>) (98) 

then 

d (P, Pi) = arccos (cos 0 cos (j) i + sin cj) sin (j) i cos (9 — 9i)). (99) 

Now choosing the north pole Pi = (0, 0,1), i.e. (j) 1 = 0 and 9\ = 0, as our fixed 
point, we get 

d(Pi,P) = 0 (100) 

and choosing the point on the equator P 2 = (1,0, 0), i.e. </> 2 = 7r/2 and 9 2 = 0, 
the distance becomes 


d (P 2 , P) = arccos (sin </> cos 9). 


( 101 ) 


To extend these two functions to R 3 , we define the Hamiltonian functions 


Hi (r (t)) 
H 2 (r (t)) 


Pi 

H 2 



= d 
= d 


Pi, 


Pi 


r(t) 

r{t) 

r(t) 

r{t) 


(j) ( 102 ) 

arccos (sin </> cos 9). (103) 


These functions can be written in terms of y/x and z/x as 


Hi ( r (0) 


H 2 (r (t)) 


Indeed, the functions 


arccos 


.\J{x/zf + {y/x) 2 (x/zf + 1, 


arccos 


, \A + (y/x ) 2 + i z /x ) 2 


hi ( x,y,z) = -, h 2 ( x,y,z) 

X 


z 

X 


are the fundamental conserved quantities of Eq. m- 


(104) 


(105) 


5.2 A Linear Poisson system 

Consider the gradient system defined by the vector field 

v ( x ) = (yz,xz,xy) (106) 

which is an unphysical version of the Euler top. The potential surfaces are 

F (x, y, z) = xyz = c (107) 
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As mentioned in [38], the geodesic flows on potential surfaces for c ^ 0 are not 
integrable, and hence one may not expect to find simple Hamiltonian functions 
for this system. However, letting c = 0 yields 

xyz = 0 (108) 

which is nothing but the non-smooth union of coordinate planes. On this surface 
consider the closed subset 


£>0, y > 0, z = 0 


(109) 


which is the first quadrant of xy— plane. Then, the Frenet-Serret frame becomes 


(t, n, b) = (k, i, j) 


( 110 ) 


Choosing Pi = (0,0) and P 2 = (1,0) on the xy —plane it is possible to define 
two distance functions 


d {Pi , P) = sjx 1 + y 2 
d{P 2l P) = y/(x - l) 2 + y 2 

whose gradients are 


Vd(Pi,P) 

Vd(P 2 ,P) 



Using the immediate condition 

j t (x(t)y(t)z(t)) = o 

we obtain the time evolution of z coordinate 

dz (t) z {t ) dx (t) z ( t) dy (t) 

dt x (t) dt y (t) dt 


( 111 ) 

( 112 ) 

(113) 

(114) 

(115) 

(116) 


That means, we can extend to the gradient vector fields in Eas. (11131) and (11151) 
to the vector fields 


ui (x) = 


u 2 (x) = 


1 


s/x 2 + y 2 
1 


(x,y,-2z) 


i - [x - l,y,-2z + 

s/d - 1) 2 + ? v 


) 


(117) 

(118) 
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which form a frame for the tangent space of the potential surface xyz = 0. Using 
the decomposition 


u 2 (x) = 


1 


(x,y,-2z) + 


1 


\J{x- !) 2 + V 2 xyj(x-l f + y 2 

it is possible find another frame consisting of gradient vector fields 

V-ffi = (x,y,-2z) 

WH 2 = (~x,Q,z) 


(~x,0,z) (119) 


for the functions 


Hi (x, y, z) = - (x 2 + y 1 - 2z 2 ) 

H 2 (x,y,z ) = ^(z 2 -x 2 ) 


( 120 ) 

( 121 ) 


( 122 ) 

(123) 


which are Hamiltonian functions of the gradient system defined by Eq. ljlOfij) . 

5.3 The Aristotelian Model of Three Body Motion 

In j3S], we showed that the Aristotelian model of three-body motion is defined 
by the gradient vector field 


v(x) = (-^- + 


:) 


x — y x — z y — z y — x z — x z — y 
where the potential surfaces are 

F (x) = a In (y — z) + b In (x — z) + c In (x — y) = K. 
We will choose the level surface defined by the constant 

K =a+6+c=l 

and parametrized by 

x — z 1 

u = - v = y- z 

y- z 2 

so that the potential surface in new coordinates will be 


In v ( u + — 


We can solve the variable v 

v = f{u) = 


»-ji =i. 


(u+i) (u~\y 


(124) 

(125) 

(126) 

(127) 

(128) 

(129) 
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and obtain the relation 


/'W 

/(«) 


u+\ u- \ 


With the parameters ( u,z ), the potential surface becomes 


X (u, z) = f ( u ) 1, Oj + z (1,1,1) 

which is a ruled surface. To obtain the orthogonal parametrization, let 
e i = —E: (1) l)l)i e 2 = —f= (2, -1, -1)) e 3 = —7= (0,1, -1) 




y/E 


V2 


and define 


y/2uf(u) f(u) 

a M = ~7s~ e2 + ~7T 

w = Vzz + -)=/ (it) fu + |) 


Then, we have 
The fundamental forms are 
( 9ij ) = 


X (u, w) = a (u) + we i. 


I«»f o 

0 1 


(hi) = 


1 


l|VF|| 


a" (■ u ) • VF 0 

0 0 


where 


vF = _ \/3f(tt) e2 + V}{f (u) + uf (u)) 


V2p (u) 


Pin) 


e 3 


a' („) = v^<"/'W + /("» e2 + />) e3 


Note also that 


VF = 


P («) 




ei x a' ( u). 


(130) 

(131) 

(132) 

(133) 

(134) 

(135) 

(136) 

(137) 

(138) 


Since the Gaussian curvature of this ruled surface is zero, the potential surfaces 
are developable surfaces, and therefore can be mapped isometrically onto the 
plane. We may choose the Frenet-Serret frame on this surface to be the Darboux 
frame 

. , , / VF a' (u) 

t,n,b ~~ 


(139) 
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of lilies of curvature which are the coordinate curves u = c\ and w = c 2 . Since 
VF, a! ( u ) and a" (u) are vectors in the plane spanned by e 2 and e 3 , we have 

VF • {a! (u) x a' 1 (it)) = 0 

and since ei is a constant vector, it is easy to see that these lines of curvature 
are also geodesics. Therefore, the coordinates (s, n, b) are geodesic distances 
in t,n, b directions, respectively. The metric in the new coordinate system 
becomes Euclidean. This allows us to write 


(n,b) = (VH 1 (n),Vtf 2 (&)) 


(140) 


where Hi and H 2 are geodesic distance functions in directions of n and b. In 
coordinates (it, v, w) , the first Hamiltonian function is 


and it is easy to see that 


f-W 

Hi (w) = / ||ei|| dt = w 

J Wn 


w = —j= (x + y + z). 
V3 


(141) 


(142) 


For the second Hamiltonian function, we have 


H 2 (u) = J \\a'(t)\\dt = J | ( tf (t) + / {t)) 2 + i (/' (t)) 2 dt (143) 


or equivalently 


H 2 ( x,y,z) 




m \ 2 

r(t)j 


+ \df (t). 


6 Conclusion 

We developed differential calculus in Frenet-Serret frame. We extended the re¬ 
sult of m for constructing Frenet-Serret frame to all dynamical systems with 
the help of eigenvectors of curl operator and a result of Chandrasekhar and 
Kendall in SS]. Considering bi-Hamiltonian structure and Jacobi identity in 
Frenet-Serret frame associated to a dynamical system, we proved that all dy¬ 
namical systems in three dimensions possess two compatible Poisson structures. 
We also presented the relation between Hamiltonian functions and Poisson vec¬ 
tors. 

Given a gradient dynamical system, we presented the geometric parameters 
of both level surface and an arbitrary curve on it. In particular, we considered, 
on level surfaces of potential function, gradient flows of restrictions of Hamil¬ 
tonian functions and proved that it is possible to find Hamiltonian functions 
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whose gradient flows on level surface have geodesic curvature zero. This result 
led us to show that Hamiltonian functions are determined by distance functions, 
namely, geodesic lenghts from an arbitrary point to two different fixed points 
on the level surface of potential function. 

Finally, by means of transformations bringing one of the components of an 
orthogonal metric to constant, we proved that finding two Hamiltonian functions 
of a gradient system is the same as constructing geodesic coordinates of its 
potential surfaces. As examples, we worked out decoupled flow of radius vector 
of a sphere, a quadratic dynamical system possessing linear Poisson structures, 
and the Aristotelian model of three body motion. 
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